Geometric interpretation of pre-vitrification in hard sphere liquids 
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We derive a microscopic criterion for the stability of hard sphere configurations, and we show 
empirically that this criterion is marginally satisfied in the glass. This observation supports a geo- 
metric interpretation for the initial rapid rise of viscosity with packing fraction, or pre-vitrification. 
It also implies that barely stable soft modes characterize the glass structure, whose spatial extension 
is estimated. We show that both the short-term dynamics and activation processes occur mostly 
along those soft modes, and we study some implications of these observations. This article synthe- 
sizes new and previous results [C. Brito and M. Wyart, Euro. Phys. Letters, 76, 149-155, (2006) 
and C. Brito and M. Wyart, J. Stat. Mech., L08003 (2007) ] in a unified view. 

PACS numbers: find packs 

I. INTRODUCTION 

Unlike crystals, amorphous structures are poorly understood on small length scales. This is apparent when one 
considers the low-temperature properties of glasses such as heat transport 1 - and the nature of the two-level systems 
leading to a linear specific heat 2 , or the statistics of force chains and stress propagation in a pile of sand 3 -. Part of the 
difficulty comes from the out-of-equilibrium nature of amorphous solids: to understand their structure and properties, 
one must also understand how they are made. This is the difficult problem of the glass or jamming transition, where 
a fluid stops flowing and rest in some meta-stable configuration. At the center of this phenomenon lies a geometrical 
question: by which processes can a dense assembly of particles rearrange, and how do these rearrangements depend 
on the particles packing? 

It is surprising that a similar question has been solved in the 70's on the apparently more complicated problem of 
polymers entanglement 4 ., where the objects considered are not simple particles but long chains forming a melt. In 
our view part of the reason for this paradox is the following: in a melt, the relaxation time scales with the length 
of the polymers. This fact can be captured experimentally and is a stringent test for theories. The situation is very 
different in glasses, where the length scales at play appear to be limited-. This fact makes it harder to distinguish 
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and compare the predictions of different theories. Nevertheless, recent numerics suggests that the length scales at 
play may not always be small. Particles interacting with a purely repulsive short-range potential display a critical 
point, corresponding to jammed packings for which the overlaps between particles vanish. Near that point, scaling 
laws characterize the microscopic structure^^, elastio^^ ' 10 ! 11 and transport properties, and relaxation in shear 
flows 13 . 

A particularly interesting observation is that soft modes, collective displacement of particles with a small restoring 
force, are abundant near this critical point^. The relation between the microscopic structure and the characteristic 
frequency and length scale of these modes was derived, and in particular the latter was shown to diverge near 
threshold^. In turn, imposing the stability of these modes led to the derivation of a non-trivial microscopic criterion 
for packing of repulsive particles^, that any mechanically stable configuration must satisfy. For infinitely fast quench 
followed by adiabatic decompression, it was observed that this criterion is marginally satisfied&i^: configurations 
generated by such a protocol are barely stable. This supported that at least for an infinitely fast quench, the 
realization of this microscopic criterion affects the dynamics, and suggested that soft modes may play a role in the 
structural rearrangements of particles. To show that this is the case in the empirically relevant situation of a slow 
quench, one would have to study a super-cooled liquid at finite temperature, and analyze soft modes, microscopic 
structure and relaxation together. This is what we perform here, using hard spheres, where interactions arc purely 
entropic and where a critical point also turns out to be present, allowing a scaling analysis. 

The paper is organized as follows. We start by illustrating the key results on the soft modes and the stability of 
packing of elastic particles using a simple model, the square lattice. In Section III, after defining the coordination 
of a hard sphere configuration, we establish a mapping between the free energy of a hard sphere system and the 
energy of an elastic network. This enables to apply all the conceptual tools developed in elastic systems to hard 
particles, in particular we derive a microscopic criterion for the stability of hard sphere configurations. In Section 
IV we present the numerical protocol we use, both in the glass and the super-cooled liquid, to identify meta-stable 
states and characterize their structural properties. In Section V we show that in those meta-stable states the stability 
criterion is saturated: configurations visited are barely stable mechanically. We confirm this observation in Section 
VI where the short term dynamics is studied. The marginal stability of the glass implies in particular an anomalous 
scaling for the mean square displacement near maximum packing, which we check numerically. In Section VII it is 
shown that only a small fraction of the degrees of freedom of the system participate to activation events where new 
meta-stable states are visited. Those degrees of freedom are precisely the soft modes present in the glass structure. 
Finally we argue that these observations support a geometric interpretation for pre-vitrification, which is presented 
in Section VIII. 
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II. A CRITERION FOR THE MECHANICAL STABILITY OF ELASTIC NETWORKS 

Studying engineering structures, Maxwell^ established a necessary criterion for the mechanical stability of elastic 
networks. The key microscopic parameter is the coordination z, the average number of interactions per particle. For 
an elastic network of springs, his criterion reads z > z c = 2d, where d is the spatial dimension of the system. The 
demonstration goes as follows. Consider a set of N points interacting with N c springs, at rest, of stiffness k. The 
expansion of the energy is: 

5E = J2 \[( 5 ^ - SRj) ■ ^l 2 + °( Sr2 )> (!) 

(ij) 

where the sum is made over all springs, n,j is the unit vector going from i to j, and 5Ri is the displacement of particles 
i. A system is floppy, i.e. not mechanically stable, if it can be deformed without energy cost, that is if there is a 
displacement field for which SE = 0, or equivalently (SRi — 5Rj) ■ n tJ = Vij. If the spatial dimension is d, this 
linear system has Nd degrees of freedom (ignoring the d(d+ l)/2 rigid motions of the entire system) and N c = Nz/2 
equations, and therefore there are always non-trivial solutions if Nd > N c , that is if z < 2d = z c . Finite stiffness 
therefore requires: 

z > 2d. (2) 

Under compression, the criterion of rigidity becomes more demanding. Here we illustrate this result in a simple 
model, but the different scaling we obtain have broader applications and are valid in particular for random assemblies 
of elastic particle o 14 i 15 . Consider a square lattice of springs of rest length a. It marginally satisfies the Maxwell 
criterion, since z = 4 and d — 2. We add randomly a density Sz of springs at rest connecting second neighbors, 
represented as dotted lines in Fig|T]), such that the coordination is z = z c + 8z. Springs are added in a rather 
homogeneous manner, so that there are not large regions without dotted springs. The typical distance between two 
dotted springs in a given row or column is then: 

I* - a/Sz (3) 

How much pressure p can this system sustain before collapsing? To be mechanically stable, all collective displace- 
ments must have a positive energetic cost. It turns out that the first modes to collapse as the pressure is increased 
are of the type of the longitudinal modes of wavelength /* of individual segments of springs contained between two 
dotted diagonal springs as represented with arrows in Fig([T]). These modes have a displacement field of the form 
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SRi = 2X sm(nia / 1*) / J I* / ae x , where i labels the particles along a segment and runs between and I* /a, e x is the 
unit vector in the direction of the line, and X is the amplitude of the mode, X = 1 for a normalized mode. In the 
absence of pressure p, springs carry no force. The energy of such a mode comes only from the springs of the segment 
and from Eq.([T]) follows SE ~ kX 2 a 2 /I* 2 . Note that these modes have a characteristic frequency: 

V m V m 

When p > 0, each spring now carries a force of order / ~ |w . The energy expansion then contains other terms 
not indicated in Eq.([T)) 15 i 17 , whose effect can be estimated quantitatively as follows. When particles are displaced 
along a longitudinal mode such as the one represented by arrows in Fig.(JTJ, the force of each spring directly connected 
and transverse to the segment considered, represented in dashed line in Fig{T]), now produces a work equal to / 
times the elongation of the spring. This elongation is simply SB 2 /a following Pythagoras' theorem. Summing on 
all the springs transverse to the segment leads to a work of order fX 2 /a. This gives finally for the energy of the 
mode SE ~ kX 2 a 2 jl* 2 — fX 2 /a, where numerical pre-factors are omitted. Stability requires SE > 0, implying that 
kcr 3 /l* 2 > /, or: 

Sz > Aif/ka) 1 ' 2 ~ (5) 

where A is a numerical constant and e is the typical strain in the contacts. This result signifies that pressure has a 
destabilizing effect, which needs to be counterbalanced by the creation of more contacts to maintain elastic stability. 
Note that Eqs. (|3l4l5[) are more general that the simple square lattice model considered here, they apply to elastic 
network or assemblies of elastic particle s 14 ' 15 as long as spatial fluctuations in coordination are limited. 

III. AN ANALOGY BETWEEN HARD SPHERE GLASSES AND ELASTIC NETWORKS 

These results on the stability of elastic networks apply to hard spheres systems. In order to see that, we recall the 
analogy between the free energy of a hard sphere glass and the energy of an athermal network of logarithm springs 18 . 
Consider the dynamics (brownian or newtonian) of hard spheres in a super-cooled liquid or glass state, such that the 
collision time among neighbor particles r c is much smaller than r, the time scale on which the structure rearranges. 
On intermediate time scale t\ such that r c << t\ << r one can define a contact network, by considering all the pairs 
of particles colliding with each other, those are said to be "in contact" (examples of contact network are shown in the 
next section). This enables to define a coordination number z. Once the contact network is defined in a meta-stable 
state, all configurations for which particles in contact do not interpenetrate are equiprobable, those configurations 



FIG. 1: Square lattice of springs with a density per particle Sz of additional diagonal springs, represented in dotted lines. 
I* ~ a/5z is the typical dimensionless distance of the segments contained between two diagonal springs on a given row or 
column. The arrows represent the longitudinal mode of wavelength ~ I* of such a segment: SRi ~ sm{iria /l*)e x , following 
the notation introduced in the text. The dashed line exemplifies the deformation of a spring transverse and directly connected 
to the segment considered, it is elongated by the longitudinal vibration of this segment. When the pressure is positive and 
contacts are under compression, this elongation lowers the energy contained in those springs. This leads to an elastic instability 
when 8z becomes smaller than a quantity proportional to the square root of the contact strain, of order f/ka. 

satisfy Yluj) — ^j'll — cr ) = 1) where O is the Heaviside function, the product is made on all contacts ij and a 

is the particles diameter, that defines our unit length. The isobaric partition function is then: 

Z = J dV]J J dRiJlQiWRi-R^-^expi^) (6) 

* (ij) 

In one spatial dimension (for a neckless of spheres), Eq.Q can be readily solved by changing variables and considering 
the gaps hij — Rj — Ri between particles in contact. The mapping is one to one and linear: 



IJ dRi cx JJ dhijS(J2 hij - (V - V ) ) (7) 

i ij ij 

where Vq is the volume of the system at p — oo. Eqs. (|6l7[ ) lead to: 



Z = H[ dhjexpiJ^L) (8) 



leading to the simple result p — k^T/{h). In higher dimensions, the situation is far more complicated in general, 
because the mapping between positions and gaps in not one-to-one, and not linear. There is nevertheless an exception 
to that rule. As was shown by several author a 20 ' 21 ! 22 , as the pressure diverges near maximum packing the system 
becomes isostatic z — > z c , see footnote^ for a sketch of the argument. As noted ir>i£, this implies precisely that the 
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number of contact is equal to the number of degrees of freedom, and that the mapping of particle positions toward 
the gaps is one-to-one. Near maximum pressure this mapping is also linear as (dRi — dRj) ■ Hij = dhij + 0(SR 2 ). One 
gets: 

Yl dR\ cx Yl dhijdiY, fijK - P(V - V )) (9) 

i ij ij 

where /y is the force in the contact ij. The volume constraint fv^ij ~ p{V ~ Vo)) generalizes for d > 1 the 

constraint 6(^2^ hij — (V — Vq)). This relation between gaps and volume can be derived as follows. In a mcta-stablc 
state, forces must be balance on all particles. As a consequence, the virtual force theorem implies that the work of any 
displacement is zero: dW = J^ij fijdhij — pdV = 0. Integrating this relation leads to the relation above . Eqs (l6l9|) 
lead to: 

Z = TT f dfc tf e--fa h «/*» r . (10) 

(y> A y >o 

and: 

f fin 

The force is inversely proportional to the average gap between particles, as we shall confirm numerically in the next 
section. The stiffness in the contact ij is then: 

kij = {k b T)/{hij} 2 (12) 

From Eq. pO)) one obtains for the Gibbs free energy Q: 

g = -hTj^Mihii)) = -hTj^HrtJ - <t) (13) 

where is the average distance between particle i and j: — (\\Ri — Rj\\)- Thus the Gibbs free energy of a hard 
sphere system is equivalent to the energy of a network of logarithmic springs. As for an elastic network, one can define 
a dynamical matrix M. by differentiating Eq. (Tig)) . M. describe the linear response of the average displacement of the 
particles to any applied force field. The eigenvectors of M. define the normal modes of the free energy^. 

When z > z c , as is the case in the glass phase (see below), Eqs. (|10tir5|) are not exact. Nevertheless, the relative 
deviations to Eq. pip can be estimated^, and are of order Sz — z — z c . Numerically these corrections turn out to be 



7 



small (smaller than 5% throughout the glass phased), and we shall neglect them. We will check this approximation 
further when we study the microscopic dynamics, see Section |VT1 Then, together with Eq.([5]) and Eq. (|12p . the present 
analogy leads to the prediction that minima of the free energy in hard spheres system must satisfy: 



where (h) and (/) are the typical gaps and forces between particles in contact, and where Eq. fTTj) was used to relate 
these two quantities. 



To study if the meta-stable states visited in the super-cooled liquid and the glass live close to the bound of Eq. (fl"4]) . 
and if the proximity of this bound affect the dynamics, we simulate hard discs with Newtonian dynamics: we use 
an event-driven code^i, particles are in free flight until they collide elastically. We use two-dimensional bidisperse 
systems of N = 64,256 and N — 1024 particles. Half of the particles have a diameter a\, which defines our unit 
length. Other particles have diameter 02 — \Act\. All particles have a mass m, our unit mass. Since for hard particles 
kbT is only re-scaling time and energy, we chose kbT as our unit of energy. Our unit of time is then <j\ \J m/k(,T. All 
data below are presented in dimensionless quantities. 

We seek to study both the glass and the super-cooled liquid phase. To generate configurations with large packing 
fractions in the glass we use the jammed configurations of£ with packing fraction distributed around 4> c ~ 0.83. At 
(j) c the particles are in permanent contact. By reducing the particles diameters by a relative amount e, we obtain 
configurations of packing fraction <fi = 4> C {1 — e) 2 . We then assign a random velocity to every particle and launch 
an event-driven simulation. This procedure enables to study the aging dynamics of highly dense systems. For 
4> < 4>q « 0.79, the system is a super-cooled liquid, and can be equilibrated. 



Computing numerically the contact network requires time-averaging on some scale t\ such that r c << t\ << t, 
where r is the a-relaxation time of the system, which we define as the time for which the self scattering function decays 
by 70%. In the super-cooled liquid a natural way to proceed would be to compute r, and chose t\ << t. Nevertheless 
this procedure is not appropriate for the aging dynamics in the glass phase, where r is not well defined, and where the 
dynamics depends on the waiting time. As an alternative protocol, we consider the self-density correlation function 




(14) 



IV. NUMERICAL PROTOCOL 



A. 



Numerical definition of meta-stable states 
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(a) (b) 
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FIG. 2: Examples of C(q,t,t u ) as defined in Eq. (|f 5 p vs t (a) in the glass and (b) in the super-cooled liquid, 
not averaged in time: 

C(q,t,t w ) = (e < *(^< t + t -)- a *< t »») i , (15) 

where the average is made on every particle j but not on time, Rj(t) is the position of particle j at time t and q is 
some wave vector. In what follows ||q|| = 2tt jo\. For all the system sizes we consider in the glass phase, and for small 
systems (for TV = 64 and to a lower extent for TV = 256) near the glass transition, we observe that C(q, t, t w ) displays 
long and well-defined plateaus interrupted by sudden jumps, as exemplified in Fig.(f5]). In the super-cooled liquid, 
those jumps are of order one, indicating that the life time of the plateaus are of order r (a few jumps de-correlate 
the structure). The existence of plateaus interrupted by sharp transitions indicates that the dynamics is intermittent, 
as previously observe d 25 ' 26 . In real space, the plateaus of C(q,t,t w ) correspond to quiet periods where particles 
are rapidly rattling around their average position. The jumps indicate rapid and collective rearrangements of the 
particles. In what follows we call "meta-stable states" those quiet periods of the dynamics. Average quantities are 
then computed in a given meta-stable state by choosing a time interval [t, t + 1\] for which the system lies in the same 
meta-stable state. We find that average quantities do not vary significantly with the location and the length of the 
time interval as long as <i >> t c . This robustness is proven for vibrational modes in particular in Annex 2. In what 
follows we chose t\ ~ 200t c . 

This protocol has the advantage to be applicable both to aging and equilibrated dynamics. On the other hand, it is 
limited to rather small systems in the liquid phase. Clearly for an infinite system C(q, t, t w ) is spatially self-averaging, 
and smooth. Already for TV = 1024 near the glass transition plateaus are hardly detectable, and our protocol does 
not apply (although it does in the glass). For such system sizes, the more traditional method (computing r from the 
decay of the smooth self-scattering function and considering some time scale t\ « r) should be used. 
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Contact force network 



A straightforward quantity to define in a meta-stable state is the average position of the particles: 



1 



t+ti 



R? = T I Ri(t')dt'. (16) 
ti Jt 

Central to our analysis is the definition of a contact force networ k 18 ' 27 ! 28 . Two particles are said to be in contact 
if they collide with each other during the time interval t\. This enables to define an average coordination number z 
as z = 2N C /N, where N c is the total number of contacts among all particles of the system. The contact force 
between these particles is then defined as average momentum they exchange per unit of time: 



7i=n co i[t 1 ] 



n=l 

where the sum is made on the total number of collisions n co /[ii] between i and j that took place in the time interval 
ti, and AP n is the momentum exchanged at the nth chock. Fig. ([3]) shows a contact force network obtained using this 
procedure, and Fig.(j4j) shows the amplitude of the contact force as a function of the average gap between the particles 
in contact. 

We define the average contact force of the network (/) as: 

1 N 

(/> = nY,wm- ( 18 ) 

Near maximum packing (/) scales as the pressure p and as the inverse of the average gap h = (hij), as implied by 
Eq. (fTTj) . The densest packing fraction O ~ 0.79 we can equilibrate corresponds to (/) « 18. For larger values of (/), 
the system is a glass. 

Note that close to maximum packing, at very large pressure, a few percents of the particles do not contribute 
to the rigidity of the structure. These "rattlers" appears in Fig.© as particles which do not exchange forces with 
any neighbors. In our analysis below we systematically remove such particles, and the procedure we use to do so is 
presented in Annex 1. 
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FIG. 3: Contact forces for N = 256, (/) = 6740 and ti = 10 4 time steps. Points represent particles centers. Contact forces 
are sketched by line segments which link particles that are in contact. The width of these segments is proportional to the force 
amplitude. This figure has already been published in the reference "On the rigidity of hard sphere glass near random close 
packing", Europhyscis Letters, v 76, 149-155 (2006) by C. Brito and M. Wyart and it is reproduced here under permission of 
the Institute of Physics Publishing (IOP). 




FIG. 4: Log-log plot of the amplitude of the contact force vs the gap between the particles for different values of (/) in a 
system with N = 256 particles. Each point corresponds to a pair of numbers (fij, (hij)) that characterizes a pair of particles 
in contact. The slashed line is a fit of the theoretical relation predicted in Eq. (|ll[) . This figure has already been published in 
the reference "On the rigidity of hard sphere glass near random close packing", Europhyscis Letters, v 76, 149-155 (2006) by 
C. Brito and M. Wyart and it is reproduced here under permission of the Institute of Physics Publishing (IOP). 



C. Normal modes of the free energy 

As shown in Eq. p3[) . the free energy in a meta-stable state can be written in terms of the average particle positions. 
It follows that it can be expanded for small average displacements. For discs (d = 2) this reads: 

«^L(0 ^ + (19) 

E of/TV" - ^)-^] 2 + o(SR 2 ) 
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where is the unit vector orthogonal to ntj . Eq. (|19[) can be written in matrix form: 



Sg = (SR\M\SR) (20) 

where |<5R) is the dN-dimensional vector 8Ri..8Rn and (tfR'^JR' 3 ) = Y£ ' ^-Rj' 3 - -A4 is the dynamical (or 
stiffness) matrix. For completeness, note that for discs it can be written as a N x N matrix whose elements Mij are 
tensors of rank d, for d = 2 this reads: 

(0 

where = 1 when particles i and j are in contact and where the second sum is made on all the particles (/) in 
contacts with the particle i. ® is the tensor product. M. describes the linear response of the average displacement of 
the particles to an external force. The eingenvectors of M. are the normal modes of the system and the frequencies are 
the square roots of these eigenvalues^,. The distribution of these frequencies is the density of states D(u>). We shall 
denote |JR W ) the displacement field of a normal mode of frequency to. These modes form a complete orthonormal 
basis {|<5R")}. 



V. MARGINAL STABILITY OF THE MICROSCOPIC STRUCTURE 

According to Eq. (|14p . minima of the free energy must have a sufficiently coordinated contact network. One may 
ask if the meta-stable states generated dynamically satisfy this bound easily, or marginally 1 —. In order to test this 
question, we prepare systems at various pressures and identify meta-stable states. Deep in the glass phase, starting 
from some initial condition typically 3 or 4 states are visited during aging on the time scales we explore. During 
aging the pressure can drop by several orders of magnitudes (indicating the possibility to obtain denser jammed 
configurations under re-compression, i.e. larger <f> c ). For each meta-stable state visited, we measure the coordination 
of the contact force network, and the average force (/). The corresponding data is presented in Fig.©, together with 
measures of the coordination in the super-cooled liquid where equilibrium is reached. 

In Fig© it appears that the fit corresponding to the saturation of the bound of Eq. (fT4)) . Sz = A{f)~ x l 2 , captures 
well all our data-points. This observation supports that the meta-stable states we generate lie close to marginal 
stability: on the time scales that can be probed numerically, the configurations visited by the dynamics have just 
nearly enough contacts to counter-balance the destabilizing effect induced by the contact forces. The situation is 
very different from a mono-disperse hexagonal crystal, for which Sz = 2 as (/) — > oo. Thus Fig© supports that, 
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FIG. 5: Log-log plot of Sz vs (/) for N = 256 or TV = 1024 particles. Each circle correspond to one meta-stable state in the 
glass phase, whereas diamonds correspond to averaged quantity among 11 meta-stable states in the super-cooled liquid. The 
slashed line is the best fit of the form Sz = A(f)~ 1 ^ 2 . This figure has already been published in the reference "On the rigidity 
of hard sphere glass near random close packing", Europhyscis Letters, v 76, 149-155 (2006) by C. Brito and M. Wyart and it 
is reproduced here under permission of the Institute of Physics Publishing (IOP) . 

at least for hard particles, there exists a fundamental difference in the mechanical stability of a glass and a crystal. 
In what follows we provide further evidences that meta-stable states lie close to marginal stability, and study some 
consequences of this property on the dynamics. 

VI. MICROSCOPIC DYNAMICS 

If a configuration is marginally rigid, then by definition it must display modes which are barely stable. In this 
section we investigate the existence of such soft modes in the free energy expansion around meta-stable states. After 
observing that these soft modes are indeed present, we show that they lead to anomalously large and slow density 
fluctuations on time scales where the system is still confined in one meta-stable state, which we refer to as " microscopic 
dynamics" . 

A. Density of States 

We compute the density of states D(ui) in meta-stables states for various pressure following the procedure introduced 
in section HV CI As the pressure is varied, following Eq. (TP2")) the characteristic stiffness and therefore the characteristic 
frequency change. It is therefore convenient to represent the density of states in rescaled frequencies lu' — Lu/(f) . 
Results are shown in Fig.^. Very similar results have been recently reported in simplified "mean field" hard sphere 
models 30 . Only the positive part of the spectrum is shown. Occasionally we observe one or two unstable modes, with 
a negative frequency, of very small absolute value. Those unstable directions may appear due to the approximation 



> 8z~<f>" 
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FIG. 6: Densities of states D'(ui') = (f)D(u) vs. rescaled frequency u)' = u>/(/) for different values of {/) in a system of 
N = 256 particles. Inset: D{u)/u> vs. uo for (/) = 18. 

we perform when computing the free energy. Alternatively, they may indicate the presence of saddles (and multiple 
configurations of free energy minima) or "shoulders" in the meta-stable state under study. 

From Fig.© we observe that: (i) there is an abundance of modes at low frequency. For all (/), D'(uj') increases 
rapidly from zero-frequency to reach a maximum at some frequency w* , before decaying again. In the inset of Fig.©, 
D(u>) is normalized by its Debye behavior Dd{ui) ~ w (plane waves would lead to a linear behavior of the density of 
states in two dimensions). No plateau can be detected at low frequency, we rather observe a peak in the quantity 
D{oj) / Dd{to), which appears at some frequency w BP significantly smaller than u>*. This indicates that for our system 
size we do not observe any frequency range where plane waves dominate the spectrum. This is confirmed by inspection 
of the lowest-frequency modes, which appear to be quite heterogeneous. Two examples of lowest-frequency modes 
are shown in Fig. ([7]) for two values of (/) . Those observations are consistent with the presence of barely stable soft 
modes in the spectrum. 

(ii) There exists a characteristic frequency uj* which scales with the pressure. Wc define uj* as the frequency at 
which D(uj) is maximum: D(ui*) — D max . Fig.© shows the dependence of uj* with the average force (/), where we 
observe the scaling: 

^-(/) 1/2 , (21) 

which holds well from the glass transition toward our densest packing, up to (/) = 10 for our system size. This scaling 
behaviour in the vibrational spectrum can be deduced from Eqs. (|4ll4p if marginal stability is assumed throughout 
the glass phase. 

Both observations (i) and (ii) bring further support on the marginal stability of the meta-stable states, previously 
inferred from the microscopic structure. 
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FIG. 7: Examples of two lowest-frequency modes for N = 1024 particles for {/) = 7.8 x 10 5 (left) and {/) = 330 (right). 
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FIG. 8: Characteristic frequency ui* as defined in the text vs average force (/). 



B. Microscopic dynamics and normal modes 

To study the microscopic dynamics, we project the dynamics on the normal modes and define for each frequency 

w. 



C u (t) = ( (SR(t + t w )\SR u ) . (5R(t w )\5R") ) tw , (22) 

where |<$R(t)) = |R(t)) — |R eq ) is the displacement field around the configuration corresponding to the average particles 
position, and where the average is made on all time segments [t wi t w + 1] entirely included in a meta-stable state. 

If the projections of the dynamics were made on longitudinal plane waves rather than on normal modes, C u (t) 
would simply correspond to the de-correlation of the density fluctuations at some wave vector, which can be probed 
in scattering experiments. Examples of C u (i) for some low- frequency modes are presented in Fig.© at two different 
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FIG. 9: Examples of Cui(t)/Cui(0) for low-frequency modes for two different average contact force {/}, deep in the glass phase 
(left) and in the super-cooled liquid phase (right). 

pressures, deep in the glass phase and in the super-cooled liquid. We observe damped oscillations for most of the 
spectrum. 

From C w (i), the amplitude A(u>) and the characteristic time t(lo) of the oscillations of a mode are readily extracted. 
The average square amplitude of the normal mode follows (A 2 (lo)) = C u (0). We define the relaxation time scale r(u>) 
as the time at which C^(i) has decayed by some fraction s: C u (t(uj)) = sC w (0). We have tried various definitions 
s = 0.3; 0.5; 0.9 and found similar scaling for the dependence of t(u>) with u>. In what follows we present the data 
with s = 0.9 where our statistic is more accurate. 

The dependence of these quantities with frequency are respectively shown in Fig. pH)) for three meta-stable states 
at different pressure. Two configurations are in the glass phase and one in the super-cooled liquid phase. In all cases, 
these quantities were computed for each mode of the spectrum. We observe that the modes display weakly damped 
oscillations, whose amplitude and period follow: 

<^ 2 H)~A> ( 23 ) 
r(u)~± (24) 

These results hold true even for the low-frequency part of the spectrum, although more scattering is found there 5 ^. 

As a consequence, our computation of D(w) gives a rather faithful distribution of relaxation time scales of the 
microscopic dynamics, supporting further the approximation we used to compute the free energy in Eq. (ll3|) . a priori 
strictly valid only at infinite pressure. This allows us to identify the peak apparent in the inset of Fig.® as the 
Boson Peak, which appears as a similar hump in Raman or neutron spectra in molecular liquid o 53 ' i 55 . Near the 
glass transition, this peak appears at a frequency significantly smaller than u>*, as shown by the inset of the Fig.© . 
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FIG. 10: (a) Average squared amplitude of the modes (A 2 (ui)) vs u) at various packing fractions in a system of 256 particles, 
both in the glass phase ({/} = 260 and (/) = 50) or in the super-cooled liquid ((/) = 18). Each point corresponds to one mode. 
The dashed line corresponds to the fit {A 2 (lo)} ~ 1/lu 2 . (b) Relaxation time t(ui) of each mode vs u for the same packing 
fraction. The slashed line corresponds to the relation t(lu) ~ 1/uj. 

C. Mean squared displacement 

In this section we use D(lu) to compute the mean square displacement around an equilibrium position inside a meta- 
stable state when (/) is varied. This quantity is directly related to the Debye- Waller factor accessible empirically 
with scattering experiments. 

We define SRi — Ri — R e L q , where R eq is the average position of particles i in a given meta-stable state as defined in 

-» 2 

Eq. (fT6]) . Assuming that the dynamics of different modes is independent, the fluctuations of particles positions (8R4 ) 
can be written as a sum of the fluctuation of all modes: 

(4 2 )=E(^ 2 H)(4H 2 ) (25) 

where A (w) is the average square amplitude of the amplitude of the mode to and SRi(co) is the displacement of 

-» -> 2 

particle i for the mode cu. We then average on all particles and define (8R 2 ) — 1/N ^2ASRi ) where TV is the system 

-* 2 

size. Using the modes normalization (SRi(u) )i = 1/N and applying Eqs.(|21|) lead to: 

(SR 2 ) ~ / ^Ld. > I ^du. (26) 

The inequality accounts for the modes with frequency between u> = and ui* that we have neglected. Accounting for 
those modes would not change our conclusion as long as the soft modes density grows sub-linearly at low frequency. 
As can be checked for the square lattice, D(u>) reaches a typical value 1/vfc (~ l/(/) for hard spheres) for to > u>*. 
This is more generally true for amorphous packing, as proven ir>ii. Using this fact, the last integral is dominated by 
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FIG. 11: Mean square displacement (SR 2 ) versus average contact force (/) for TV = 1024 (circles) and N = 256 (squares) 
particles. Diamonds correspond to the super-cooled liquid phase and were computed for a system with N — 256 particles. 
Slashed line corresponds to the best fit agreeing with our prediction (SR 2 ) ~ (/) -3 ^ 2 - 



the lowest bound and one gets: 

(6R)>^~(.f)- 3/2 ~h^ (27) 

which holds in any dimension d > 2 (with corrections of order h 2 logiV for d = 2 due to plane waves). We have used 
the scaling of the frequency scale oj* confirmed in Fig. (JBJ) . In crystals, the fluctuations around a particle position is of 
the order of the inter-particle gap h: (SB 2 ) ~ h 2 (with log TV corrections in two dimensions). Eci. (f27|) shows that, near 
maximum packing, the amplitude of particles motions is infinitely smaller in the crystal than in the glass. Because of 
the marginal stability of the glass, these fluctuations have an anomalous scaling with the packing fraction. 

To check numerically this prediction, we consider various meta-stables states. In each of them, we measure R^ q 
and the mean square displacement around the equilibrium position: (SR. 2 ) = (1/TV ^ i <5i?f (i)} t x , where the average 
is made on the time interval t\. Fig. (|ll[) show this quantity for various packing fraction. Our numerical result agrees 
well with our prediction (5R 2 ) ~ (/)~ 3//2 ~ h 3 / 2 throughout the glass phase. 



VII. a-RELAXATION 



One long-lasting challenge in our understanding of the glass transition is the elaboration of a spatial description 
of activated events, the rare and sudden rearrangements of particles corresponding to jumps between meta-stable 
states. These events are collective rearrangements of particles, but the cause and the nature of this collective aspect is 
unknown. Our observation that the glass structure is marginally stable suggests that the softest, barely stable modes 
may play a key role in the activated events that relax the structures. In what follows we investigate this possibility 
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by projecting the sudden rearrangements on the normal modes of the free energy. 

A. Aging 

During aging, sudden rearrangements, or "earthquakes", appear as drops in the self scattering function, see Fig.((2]- 
a). Such earthquakes correspond to collective motions of a large number of particles, and have been observed in 
various other aging systems, such as colloidal paste or laponite 3 -^, and in Lennard-Jones simulation o 26 ' 32 ! 33 . Even for 
our largest numerical box of N — 1024 particles, deep in the glass phase these events generally span the entire system. 
Examples of earthquakes in real space are shown in Fig. (|14[) . 

To analyze these displacement fields, we measure the average particle positions and the contact network in the 
meta-stable state prior to the earthquake, and compute the normal modes of the free energy. We also compute 
the earthquake displacement |^R e ) defined as a difference between the average particles position in two successive 
meta-stable states I and m: \SH e ) = |R' m ) — |R'). We then project |<5R e ) on the normal modes and compute 
c w = (<5R e |<5R w )/(<5R e |<5R e ). The cjs satisfy X^ c ^ 2 = 1 since the normal modes form a unitary basis. To study 
how the contribution of the modes depends on frequency, we define: 

<?H = (cJ) (28) 

where the average is made on a small segment of frequencies [w,v + duj]. Fig. (fT2l a) shows g(u>) for the earthquake 
shown in Fig.^a). The average contribution of the modes decreases very rapidly with increasing frequency, and most 
of the displacement projects on the excess-modes present near zero-frequency. This supports that the free energy 
barrier crossed by the system during a rearrangement lies in the direction of the softest degrees of freedom. 

To make this observation systematic, introduce the label i to rank the c's by decreasing order: c\ > C2... > c 2 n. 
We then defined a fci/2 such that: 

fcl/2 

E C NV2 (29) 

Physically, k 1 / 2 is the minimum number of modes necessary to reconstruct 50% of the displacements relative to the 
earthquake. Fig. (fT5|) shows F 1 / 2 = k 1 / 2 /(2N) for the 17 cracks studied and indicates that 0.2% < F 1 / 2 < 2% for 
all the events studied throughout the glass phase. We thus systematically observe that the extended earthquakes 
correspond to the relaxation of a small number of degrees of freedom, of the order of 1% of the modes of the system. 

In Fig. (fT4")) we illustrate the spatial consequence of our analysis. Two examples of earthquakes at different packing 
fractions are compared with the linear superposition of the 1% of the modes that contribute most to them. The 
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FIG. 12: Straight curve: D' (to 1 ) / D' max vs. ui/(f). Both D(u>) and {/} are computed in the meta-stable state prior to the 
earthquake shown in the Fig.([2]-a). Dotted curve: g(u) as defined in the text re-normalized by its maximum value g{oj') / g max 
vs. u)/(f). This figure has been originally published in the reference "Heterogeneous dynamics, marginal stability and soft 
modes in hard sphere glasses", J. Stat. Mech., L08003, (2007), by C. Brito and M. Wyart. 
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FIG. 13: Fi/2 vs. {/} for N = 256 (circles) and N = 1024 (diamonds) particles. This figure has been originally published in the 
reference "Heterogeneous dynamics, marginal stability and soft modes in hard sphere glasses", J. Stat. Mech., L08003, (2007), 
by C. Brito and M. Wyart. 



similarity is striking: the complexity of the structural relaxation is indeed contained in the soft degrees of freedom of 
the system, along which yielding occurs. Thus, in this regime only a small fraction of the degrees of freedom of the 
system participates in the relaxation of the structure. 



B. Structural Relaxation in the equilibrated super-cooled liquid 

We equilibrate the system for a range of density 0.77 < 4> < 0.786. Also in this regime, the dynamics is heterogeneous 
in space and in time and sudden rearrangements still occur on time scales of the order of t— . An example of this 
rearrangement, that can be identified as a drop in the self-scattering function, is shown in the Fig.([2]-b). In real 
space, this displacement corresponds to a collective event, as one can observe in the examples of the Fig. ([Pol - left) and 
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FIG. 14: Above: two examples of earthquakes in the glass phase for different average contact force (/) for N — 1024 particles. 
Displacements were multiplied by four for visibility. Below: projection of earthquakes on the 1% of the normal modes that 
contribute most to them. 

Fig. (JTHJ-above) . To study these events in an equilibrated super-cooled liquid, we extend the procedure used in the 
aging regime: we identify the meta-stable states visited by the dynamics and compute their averaged configuration. 
We then define the normal modes in the meta-stable state and the displacement field corresponding to the relaxation 
events. For each relaxation event, we compute F 1 / 2 - 

We start with a system with N — 64 particles. For each packing fraction, F 1 / 2 is computed for six relaxation 
events. Then this quantity is averaged on all events. The result of these (-F1/2) are shown in Fig. tfTS"]) as a function 
of the packing fraction. We find that (-F1/2) < 5% for all cf) studied, supporting that only a small fraction of the 
low-frequency modes contribute to the structural relaxation events also in this region of the super-cooled liquid. This 
fraction decays significantly as <p get closer to 4>q, suggesting a rarefaction of the number of directions along which the 
system can yield near the glass transition. Fig. (|16p exemplifies this conclusion: for this particular case, the relaxation 
event shown on the left projects almost entirely on one normal mode, shown on the right in this same figure. This 
mode turns out to be the lowest-frequency normal mode of the free energy. 

To study finite-size effects, we measure -F1/2 for twelve relaxation events at each of the five packing fraction using 
TV = 256 particles. Results are shown in Fig. ([T5|) . Finite size effects are present, and (-F1/2) appears to be roughly 0.5% 
higher in the larger system for all packing fractions. Most of this difference in behavior is explained by the observation 
that the glass transition occurs at smaller packing fraction in the N = 64 system, as previously observed^. The inset 
of Fig. (jTTJ) shows that this is the case in our system as well. If (-F1/2) is plotted as a function of relaxation time, 
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FIG. 15: (F1/2) vs <j> for two different system sizes. This figure has been originally published in the reference "Heterogeneous 
dynamics, marginal stability and soft modes in hard sphere glasses", J. Stat. Mech., L08003, (2007), by C. Brito and M. Wyart. 
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FIG. 16: Left: displacement field for a system with N = 64 particles corresponding to a relaxation event. Arrows were 
multiplied by 1.2. Right: normal mode that contains 80% of the projection of the real displacement field. This normal mode 
has the lowest frequency of the spectrum. This figure has been originally published in the reference "Heterogeneous dynamics, 
marginal stability and soft modes in hard sphere glasses", J. Stat. Mech., L08003, (2007), by C. Brito and M. Wyart. 



as in Fig. (|17p . the curves become similar for the two systems and (F 1 / 2 ) is systematically smaller for a system with 
./V = 256 particles. Thus, even for larger systems, collective rearrangements relaxing the system are "soft": they 
project mostly into a small portion of the vibrational spectrum. We verify spatially this observation in Fig. (|18p . 
Three examples of relaxation events at different packing fractions are compared with the vector field which is a linear 
superposition of the modes that contribute most to them. This relation between soft modes and relaxation has been 
recently supported by the observation that regions where structural relaxation is likely to occur, said to have a high 
"propensity" , also display an abundance of soft modes^. 
Interestingly, the soft modes that characterize marginally stable structures are in general rather extended objects, 
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FIG. 17: {-F1/2) vs t. Inset: a-relaxation time vs <j>. System sizes are indicated in the legend. This figure has been originally 
published in the relerence "Heterogeneous dynamics, marginal stability and soft modes in hard sphere glasses", J. Stat. Mech., 
L08003, (2007), by C. Brito and M. Wyart. 
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FIG. 18: Above: Relaxation events in the super-cooled liquid for different cj>, indicated in the figure, and N = 256 particles. 
Displacement fields are rescaled by 4, 1.5 and 1.2 respectively for visibility. Below: projection of the relaxation events on 
the normal modes that contribute the most. The fraction of the total number of modes used is indicated in each figure, and 
corresponds to the fraction necessary to recover 50% of the relaxation event. As indicated in Fig, (|15p . this fraction of modes 
tend to increase as tf> decreases. 



as can be observed from the examples presented here. Theoretically this is what one expects both in the square 
lattice, as justified by Eq.((3]), and in amorphous packin g 14 i 15 . In this light it does not seem surprising that activated 
events are collective. 
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VIII. A GEOMETRIC INTERPRETATION OF PRE- VITRIFICATION 

A. Pre-vitrification 

We have shown, both from its microscopic structure and microscopic dynamics, that the hard sphere glass lies close 
to marginal stability. In this section we propose an explanation for this observation. This requires a single assumption, 
namely that the viscosity increases very rapidly when meta-stable states appear in the free energy landscape. In the 
logarithmic representation of the plane coordination vs the typical gap between particles in contact (Sz, h), there 
exists a line, corresponding to the equality of Eq. (|14p , which separates a region where configurations are stable and 
unstable, as sketched in Fig. ()19|) . At any packing fraction (j>, equilibrium configurations correspond to a point in the 
(Sz, h) phase diagram. As <fi is varied equilibrium states draw a line in this plane, represented by the dashed line (red 
online) in Fig. (|19[) . At low <f>, gaps among particle are large and configurations visited are unstable. As 4> increases, 
the gaps narrow, and configurations become eventually stable. This occurs at some iponset when the equilibrium 
line crosses the marginal stability line. At larger <f>, the viscosity increases sharply, so that on our numerical time 
scales equilibrium cannot be reached deep in the regions where meta-stable states are present. As a consequence, 
the system falls out of equilibrium at some 4>o close but larger than 4> nset ■ Configurations visited must therefore lie 
close to the marginal stability line, as represented by the dotted line in Fig. (TT5]), since more stable, better-coordinated 
configurations cannot be reached dynamically. 

In this view, <f> nset corresponds to the onset temperature, where activation sets in and the dynamics becomes 
intermittent. When intermittency appears, the a-relaxation time scale r is still limited, and has increased roughly 
of one order of magnitude from the liquid state. This is consistent with the observation that the configurations we 
probed in the super-cooled liquid, for which r is larger but still limited, are already stable: the free energy expansion 
has in general a positively-defined spectrum. Note that at 4> nset, Sz ^ 0, and the characteristic length of the soft 
modes 1* is finite. We think of those modes as involving a few tens of particles. 

B. Ideal glass and random close packing 

To put our work in a broader context it is useful to think about the phase diagram of Fig|19lwith an extra dimension 
added. For any configuration one can associate the packing fraction ip c corresponding to the jammed packing that 
would be obtained after a rapid compressio n 30 ' 36 i 37 . At equilibrium (f> c is an increasing function of pressur e 38 ' 39 . 
In the three-dimensional phase diagram (h, Sz, ip c ), marginality is now represented by a surface. Its main feature, 
the scaling relation between coordination and typical gaps expressed in Eq. (|I4|) , holds irrespectively of the value of 
4> c according to our theoretical analysis, in agreement with the data presented in Fig. ([5]). Marginality and related 
properties are therefore adequately discussed in the more simple two-dimensional phase diagram presented in Fig |19l 
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FIG. 19: Phase diagram for the stability of hard sphere configurations, in the coordination 8z vs average gap h plane. The 
marginal stability line delimits stable and unstable configurations. The dashed line correspond to equilibrium configurations 
for different (f>. As <f> increases, h decreases and the two lines eventually meet. This occurs at the onset packing fraction 4> m set, 
where dynamics become activated. At larger <j>, viscosity increases sharply as configurations visited become more stable. For a 
finite quench rate the system eventually falls out of equilibrium. More stable and more coordinated regions cannot be reached 
dynamically, and as <f> is increased further, the system lives close to the marginal stability region, as indicated in the dotted 
line. The location of the out-of-equilibrium trajectory depends on the quench rate. In the limit of very rapid quench, the 
out-of-equilibrium line approaches the marginal stability line. 

Some other aspects of the dynamics and thermodynamics of hard spheres nevertheless benefit from introducing the 
extra dimension 4>c- 

Ideal glass: For molecular glasses the presence of an ideal glass transition where the configurational entropy vanishes 
at finite temperature has been proposed and is still debate d 40 ! 41 . For hard spheres this view implies that the viscosity 
diverges at some finite h > 30 i 36 , for which the equilibrium curve <f> c {h) reaches a constant value. Our work does not 
address the issue of the existence of an ideal glass, but it supports that if it exists, it is not responsible for the slow-down 
of the dynamics in the pre- vitrification region we can access empirically, neither in the aging regime we could observe 
in the glass, since such a scenario would not explain the marginal stability of the microscopic structure we observe. 
Some authors have used diverging fits of the relaxation time to argue in favor of the opposite view- 3 ^. Nevertheless, 
establishing the existence of an actual divergence from such fits is questionable even in molecular liquids 4 ^ where the 
number of decades of viscosities accessible is two to three times larger. 

Random Close Packing: Empirically it is observed that for various protocols of compression (such as pouring metallic 
balls in a container), the final packing fraction obtained is <p c ss 0.64 for mono-disperse hard spheres. This fact can be 
expressed as follows: for each protocol one can associate a line 5z(h),(j) c (h) characterizing the configurations visited 
during compression. Isostaticity implies 6z — > as h — > 0. Furthermore, for a wide class of protocols 4> c {h) — > 0.64 as 
h — ► 0. The explanation for this observation is debated 3 ^. An interesting hypothesis is that 4> c w 0.64 corresponds to 
the limit reached by infinitely rapid compressions. If typical protocols are fast in comparison with the relevant time 
scales of the dynamics, they should generate packings with a similar packing fraction. 
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IX. CONCLUSION 

We conclude by a brief summary of our results and a few remarks. We have derived a geometric criterion for 
the stability of hard sphere configurations, and we have shown that in a hard sphere glass this bound is nearly 
saturated. This supports that pre- vitrification occurs when the coordination is sufficiently large to counter-balance the 
destabilizing effect of the compression in the contacts. Nearly unstable modes are collective displacement fields, whose 
spatial extension is governed by the coordination. Once meta-stable states appear in the free energy, activation occurs 
mostly along a small fraction of these soft modes. This observation supports that these modes are the elementary 
objects to consider to describe activation. It also implies that structural relaxation must be cooperative, since the 
soft degrees of freedom are collective. 

We have observed that less and less modes participate to the structural relaxation as the packing fraction increases 
near (j) nset- It is tempting to speculate that, as the number of degrees of freedom allowing relaxation is reduced, the 
size of the cooperatively rearranging regions grows to eventually saturate at the extension of the softest modes I*. 
Nevertheless, a quantitative description of the relationship between soft modes and dynamical length scale remains 
to be built and tested. Other factors, such as the possible presence of locally favored structure of high coordination 
or some other spatial heterogeneities of the structure, may also have to be taken into account. 

Our analysis of the structural relaxation at equilibrium applies to the pre- vitrification region, corresponding to the 
intermediate viscosities that can be probed numerically. Similar time scales are accessible experimentally in shaken 
granular matter and colloidal glasses. Our work does not address the behavior of the equilibrium dynamics for very 
large packing fraction. As a consequence, it is possible that at much larger viscosities than those we probed, in 
particular near the glass transition in molecular liquids, our observations on the nature of the structural relaxation 
may not apply, and soft modes may play no role in the dynamics. Nevertheless several observations support that soft 
modes and dynamics are related even for those large viscosities. In particular, the intensity of the boson peak, which 
indicates the presence of soft modes in the spectrum, strongly correlates with the glass fragility 4 ^, a fact which is not 
captured by current theories of the glass transition. 

Finally, our geometric approach to pre-vitrification is consistent with Goldstein 4 ^ views, who proposed 40 years 
ago that the glass transition is related to the emergence of meta-stable states in the energy landscape. Other more 
recent descriptions of the glass transition, such as the mode coupling theory (MCT) of liquids 45 , make a similar 
predictio n 41 i 46 ' 47 i 48 , and it is interesting to compare this approach to ours. Here wc indicate several differences and 
analogies in the respective conclusions: (i) in MCT the predicted location of the elastic instability corresponds to 
the onset packing fraction 49 . This is consistent with our observation that when the dynamics becomes intermittent 
(</> > 4>onset)<, the configurations visited have in general a positively defined spectrum, displaying no unstable modes. 
This is also supported by previous results showing that the dynamics is dominated by activation in this parameter 
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ranged. Nevertheless, MCT predicts diverging time scales^ and dynamical length scales^ at the onset packing 
fraction, which are not observed. Fitting empirical data with such divergences^ leads to a critical packing fraction 
4>mct significantly larger than 4> nset- The interpretation of the extra fitting parameter <^mct, and its relation with 
the free energy landscape, is at present unclear, (ii) In MCT the dynamics is computed via a resummation of a 
perturbation expansion in the non-linear interaction among modes, around a point where plane waves are un-coupled. 
In our case, we use a variational argumen t 14 i 15 to capture the properties of the linear soft modes whose stability 
is at play. This argument applies as well in covalentl and attractive glasse o 7 ! 56 . This leads to an estimation of a 
length scale I* characterizing soft modes, which depends on the coordination. This length scale has not yet found 
a correspondence within MCT, where non-trivial length scales appear from the dynamics^ but diverge near the 
clastic instability, unlike I*, (iii) In our approach, the key microscopic parameters determining the location of the 
transition are coordination and pressure. In MCT, an important parameter is the area under the first peak of the 
pair correlation function^.. These two views bear similarities, as the later quantity can be considered as a rough 
measure of coordination. It remains to be seen if MCT can capture the critical behavior of the marginal stability line 
observed at very large pressure. Exploring this possibility may clarify the physical meaning of the approximations 
that characterize MCT. 
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APPENDIX A: DETERMINATION OF THE RATTLERS 

Near maximum packing, a few percents of the particles are trapped in a large "cages" on which they apply a minus- 
cule force in comparison to the typical contact forces in the system. Such particles, called rattlers, do not participate 
to the rigidity of the structure: if removed, stability is still achieved. When we compute e.g. the coordination of the 
microscopic structure, we do not take these particles into account. 

To identify rattlers we measure the average number of shocks per contact for each particle. We compute how many 
shocks n s hoc and how many contacts n c each particle has during the interval of time t\ and define: /* = n s h oc jn c 
if n c > 2 and /* = otherwise. This quantity is normalized by the average number of shocks per contact that all 
particles have during t\. F* = N s hoc/N c , where N s h oc in the total number of shocks and N c is the total number of 
contacts in the system. We then plot the distribution of f*/F* for different packing fractions, Fig. (f2T)]) . At large (/), 
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we observe the emergence of a peak near zero. When (/) is intermediate, (/) = 5.2 x 10 3 and (/) = 9.2 x 10 2 , the peak 
vanishes. This peak corresponds to the rattlers. In this work we consider that all particles for which f* JF* < 2% are 
rattlers. This criterion is represented by the arrow in the inset of the Fig. (|20|1 . To check the robustness of our results, 




FIG. 20: Histogram of the distribution of f* /F* (see definition in the text) for various average force. 



we test if the relation between the excess of coordination Sz and the average force (/) depends on this criterion. We 
vary the threshold bellow which we consider a particle as a rattler and plot in the Fig. pTj) the comparison between 
3 different criteria: f*/F* < 0.01 < 0.02 < 0.05. We observe that relation Sz = Ai(f)^ 1 / 2 holds irrespectively of 
the criterion. It fails when the rattlers are not removed of the analysis. In this case, for high values of (/) one finds 
Sz < 0. 
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FIG. 21: Sz vs (/) for 3 different criteria of definition of rattlers as explained in the text. The legend "all particles" indicates 
that no rattlers are excluded for this measure. Dotted curve: fit of the relation Sz — A(f)~ 1 / 2 . 
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APPENDIX B: PERSISTENCE OF THE NORMAL MODES IN A META-STABLE STATE 

Here we show that our numerical computation of the normal modes is robust to different choices of time intervals, 
as long as they lie in the same meta-stable state. To achieve that we compute the normal modes for two distinct, 
non-overlapping time-intervals. |<5R" ) denotes the normal modes of frequency ui, computed on some time interval 
labeled t a . We then compute the matrix of scalar product: 

C u ^ = {8B!i a \SB£). (Bl) 

If the two sets of normal modes are identical, C is the identity matrix. In general this must not be exactly true, as 
shown in Fig |22i since our protocol requires time-averaging and is therefore noisy to some extent, and also because 
some non-trivial dynamics may still occur within meta-stable states. Our observations below show that those effects 
are small, even if the two time- intervals considered are separated by a time scale of the order of the life-time of 
meta-stable states. To quantify the difference between C and the identity matrix, we follow the procedure we used 
before to compare a relaxation event to the normal modes of the structure. We define Fi/ 2 (oj) as the minimal fraction 
of normal modes computed on tb sufficient to represent 50% of a normal model of frequency uj computed on t a . We 
then define (F 1 / 2 ) as the average of F 1 / 2 (ui) on the 20 lowest-frequency modes computed on t a . (Fi/ 2 ) is 1/2N if 
C is the identity matrix, and should be small if our procedure is robust to different choice of time-interval. This is 
indeed the case: in the super-cooled liquid ((/) = 18) we find (-F1/2) = 0.4%, which is small for all practical purposes 
discussed in this article. For this measure the time intervals lasted t\ = 500r c , and the time separation between t a 
and tb was 10 4 r c , which is of the order of the relaxation time t«3x 10 4 t c . For the glass ((/) = 6700) F 1 / 2 (u)) is 
close to 1/2N, as is obvious from Fig |2"21 
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FIG. 22: Matrix of scalar product C as defined in the text for the 100 lowest-frequency modes for N = 256. The colorbar 
indicates the value of the scalar product. Time intervals lasted t\ = 500r c . Left: {/} = 6700. The intervals t a and tb are 
separeted by 10 5 r c . Right: (f> = 0.782. The intervals t a and ti are separeted by 10 4 r c , which is of the order of the relaxation 
time. 
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